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Using new four-loop results for the heavy quark vacuum polarization and new data for bottom 
quark production in electron-positron annihilation, an update on the determination of charm- and 
bottom-quark masses through sum rules has been performed. The previous result for the charm- 
quark mass, m c (3 GeV) = 0.986(13) GeV, based on the lowest moment, is supported by the new 
results from higher moments which lead to consistent values with comparable errors. The new value 
for the bottom quark, m(,(10 GeV) = 3.610(16) GeV, corresponding to mb(mb) = 4.163(16) GeV, 
makes use both of the new data and the new perturbative results and is consistent with the earlier 
determination. 
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I. INTRODUCTION 

The precise determination of charm and bottom quark 
masses has always been an important task both for the- 
ory and experiment. The most precise values have been 
obtained [l[ from an analysis of the ITEP sum rules Q 
(for reviews see Refs. B 0, 1), combining data for the 
heavy quark production cross section in electron-positron 
collision with dispersion relations and a four-loop evalu- 
ation of the vacuum polarization induced by the heavy 
quark current. In this letter, we present an update of 
these results. We will include data recently published 
by the BABAR collaboration Q and make use of new 
perturbative results which replace the estimates for the 
four-loop coefficients of higher moments used in the ear- 
lier publications. 



II. ANALYTIC RESULTS 

Our determination of the heavy quark masses follows 
closely Refs. 0, 0, S|- It is based on the direct compar- 
ison of the theoretical and experimental evaluations of 
the contributions to the derivatives of the polarization 
function LIq^ 2 ), the former evaluated in perturbative 
QCD, the latter through moments of the measured cross 
section for heavy quark production in electron-positron 
annihilation. Using dispersion relations, the moments of 

Mn^J^RQis), (1) 

can be related to the derivatives of the vacuum polariza- 
tion function at q 2 = 



Mn 



12n 2 



(2) 



In its domain of analyticity LlQ(g 2 ) can be cast into the 
form 



n Q (g 2 ) = 



'16tT 2 



(3) 



n>0 



with z = q 2 /{Am'o). Here wlq = rriQ(fi) is the heavy 
quark mass with charge Qq in the MS scheme at the 
scale fj,. The coefficients C n depend on a s and on 
the heavy quark mass through logarithms of the form 
l mQ = ln(TOg(^)/^ 2 ). Equating theoretically calculated 
and experimentally measured moments, the heavy quark 
mass is given by 



1 (*? Q Cn\ 1/(Pn) 
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As a perturbative series the coefficients C n can be written 

as 

'a s (n)\ 3 f- (30) ^(31), ,/7(32)i2 

\ n ' L ra Q i <~'n <- mQ 



(5) 



C^ll 



q 2 =0 



The terms of order a s were evaluated up to n = 8 in 
Refs. @, [13, El (and recently in Refs. (P3 . [l3l | even up 
to n=30). The four loop contributions to Co and C\ 
were calculated in Refs. [lj, [lij]. For the higher moments 
the analysis of [H was based on estimates for Cn°^ with 
n = 2, 3, 4, which lead to an additional uncertainty in 
the mass determination. Recently the exact results for 
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n 


1 


2 


3 


4 


charm 


-5.6404 


-3.4937 


-2.8395 


-3.349(11) 


lower | upper limits 




-6.0 || 7.0 


-6.0 || 5.2 


-6.0 || 3.1 


bottom 


-7.7624 


-2.6438 


-1.1745 


-1.386(10) 


lower 1 1 upper limits 




-8.0 || 9.5 


-8.0 || 8.3 


-8.0 || 7.4 



TABLE I: New results for the coefficients cL 30) m comparison 
with previous upper and lower limits as used in Ref. For 
less precise numerical results of C n for n — 3 and n = 4 see 
Ref. 0. 



the second |16| and third [17| moments were obtained. 
Combining these coefficients with additional information 
on the threshold and the high-energy behaviour and us- 
ing the analyticity of TlQ(q 2 ) and Pade approximations, 
fairly precise numerical results were obtained [HI for the 
higher coefficients up to n = 10. (For an earlier analysis 
along similar lines see Ref. [HI].) For the lowest four mo- 
ments the four loop coefficients Cn°^ are listed in Tab. U 
both for the charm and the bottom quark. All other co- 
efficients relevant for n — 1 to 4 can be found in Tabs. 4 
and 9 of It should be emphasized that these results 
are well within the estimates used in the analysis of [![, 
also shown in Tab. fl] The impact of these new results on 
the quark mass determination will be studied below. 



III. BOTTOM PRODUCTION CLOSE TO 
THRESHOLD 



The determination of the bottom quark mass, as per- 
formed in 0, 0] relies heavily on the precise measure- 
ment of R = <r(e + e~ — > hadrons) /(Tpt (with cr pt = g" ), 
which enters the moments as defined above. Specifi- 
cally, it is the contribution from the heavy quark cur- 
rent denoted as Rb with the light quark contribution sub- 
tracted. It is convenient to split the integration region 
into three pieces: The lowest region covering the nar- 
row resonances, an intermediate "threshold" region be- 
tween 10.62 GeV and 11.24 GeV, and the perturbative 
region above 11.24 GeV, where the measurement is re- 
placed by the perturbative QCD prediction. The choice 
of 11.24 GeV corresponds to the upper end of the en- 
ergy range covered by a CLEO measurement more than 
20 years ago [2(j. It also coincides approximately with 
the energy reach of a recent BABAR measurement [f|. 
In the analysis of [H, T(45) with its mass Afx(4S) = 
10.5794(12) GeV and width r T(4S ) = 20.5 MeV has been 
considered together with the three lower, narrow reso- 
nances and thus the continuum part of the bottom-cross 
section was taken from 10.62 GeV upwards. Until re- 
cently the only measurement in the threshold region has 
been the one from the CLEO collaboration, which quotes 
a systematic error of about 6%. No radiative corrections 
had been applied. In Ref. [l[ it has been argued, that 



a normalization factor 1/1.28 is necessary to reconcile 
these data with more recent and more precise CLEO re- 
sults below the T(4S')-resonance and with perturbative 
QCD at the high end. These "rescaled" data were the 
basis of the subsequent extraction of the bottom quark 
mass. However, in view of these uncertainties an overall 
systematic error of 10% was attributed to the contribu- 
tion of the moments from this region. Thus, although 
this contribution to the moments is relatively small, its 
impact on the error was larger or equal than the one from 
the other two regions combined. 

Recently a measurement of Rb in the energy region be- 
tween 10.54 GeV and 11.20 GeV was performed by the 
BABAR collaboration with significantly improved statis- 
tics and with a correlated systematic error between 2.5% 
and 3% Q. In principle this should allow an indepen- 
dent determination of the contribution to the moments 
with significantly reduced systematic error. However, no 
radiative corrections were applied to the published data 
and the radiative tails of the four lower T resonances were 
included in the quantity denoted Rb- In the following we 
describe the procedure used to obtain the contribution 
to the moments from these data. 

In a first step we subtract the radiative tail of the 
T(15), T(25) and T(3S) resonances, which is explicitly 
given in Ref. jf|. Subsequently we subtract the radiative 
tail of the T(45) resonance. For the resonance shape 
we use a Breit-Wigner function with an electronic width 
of r ee (T(4S 1 )) = 0.272 keV and an energy- independent 
total width of r tot (T(45)) = 20.5 MeV |lj|. For the ra- 
diator function G(z) we take the functional dependence 
as used in [22|, based on the resummed NNLO result of 



G{z) = (3(1 - zf~ x e 5 ^ F (6% +s 
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The remainder a corresponds to the continuum cross 
section distorted by initial-state radiation (ISR) and 
modified by vacuum polarization. It is related to <7, the 
cross section without ISR, through 



a(s) 



dz G(z) a(sz) . 



(8) 



where the lower bound of the integration is given by 
zq = (10.62GeV) 2 /s corresponding to the point where 
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FIG. 1: Comparison of rescaled CLEO data for Rb with 
BABAR data before and after deconvolution. The black bar 
on the right corresponds to the theory prediction [24} . 



the continuum cross section (after subtraction of the 
T(4s) resonance) vanishes. 

Given a, we can solve for a in an iterative way as fol- 
lows: Let us define SG(z) = G(z) — 5(1 — z) and evaluate 
a successive series of approximations, 



dz5G(z)<Ji-i(sz) 



(9) 



using as starting point cr = <j. The difference between 
<Ji and tr can be estimated by evaluating Eq. ([5]) with 
in place of a. After five iterations the resulting function 
differs from a by less than 0.5%. 

Finally, the effect of the vacuum polarization must be 
taken into account and the result is normalised relative 
to the point cross section, 



Rb 



3s 



47ra 2 (s) 



(10) 



The integration region in Eq. ([5]) covers the energy 
interval between 10.62 GeV and 11.24 GeV, whence a 
constant value (a/a(s)) 2 — 0.93 has been adopted. 

In Fig. [T]we show the BABAR data Q (after subtrac- 
tion of the radiative tails of T(15) to T(4S)), together 
with Rb after deconvolution of ISR and correcting for the 
running of a(s). We also show the CLEO data [2fJ after 
the aforementioned rescaling. 

It is now straightforward to evaluate the contribution 
to the moments. The result is listed in Tab. lU and 
com par ed to our earlier analysis based on the CLEO re- 
sult [20| . The error of this new result is dominated by the 
correlated systematic error of the BABAR measurement 
which amounts to about 3.5%. In addition we use a 2% 
error for the uncertainty from the details of the match- 
ing between the tail of T(AS) and the continuum around 
\/s = 10.62 GeV, which we add in quadrature. 



n 


12 3 4 


M dat x 10O+ 1 ) 

n, old 

A/fdat x in(2n+l) 

■ /Vl n, new ^ ±XJ 


0.296(32) 0.249(27) 0.209(22) 0.175(19) 
0.287(12) 0.240(10) 0.200(8) 0.168(7) 


W^ncw X 10^+1) 


4.592(31) 2.872(28) 2.362(26) 2.170(26) 



TABLE II: Moments in (GeV)" 2n for the bottom quark sys- 
tem from the threshold region (10.62 GeV, 11.24 GeV) from 
Ref. [lj] (old) and this letter (new). Also the new total exper- 
imental moments are given. 



As is evident from Tab. HIl the agreement between old 
and new result is remarkable giving additional confidence 
in the procedure used in Ref. [1J. The new experimental 
input and the new information on the coefficients 
lead to a significant reduction of the error on mj, as 
shown below. 



IV. QUARK MASSES 

In the absence of new data the analysis of m c will be 
based on the moments listed in Tab. 6 of Ref. [l[. As 
emphasized earlier P, 0] it is convenient to consider as 
primary quantity the running quark mass at scale 3 GeV. 
This is the natural scale for the sum rule (corresponding 
roughly to the charm threshold) and, as a consequence 
of the smaller strong coupling constant, the perturbative 
series exhibits a more stable behaviour. 

If not stated otherwise, all input parameters and as- 
sumptions are identical to those of Ref. [l|. In partic- 
ular we adopt a s (Mz) = 0.1189. The new results and 
the corresponding errors are listed in Tab. IIIIl Com- 
pared to [l|] , the shift induced by the analytic results for 
Cn°^ amounts to 3, 4 and 8 MeV for n=2, 3 and 4 re- 
spectively. The results for all four moments are nicely 
consistent, and the three lowest moments exhibit compa- 
rable errors, ranging between 13 MeV and 17 MeV. Note, 
that the relative composition of the experimental input 
varies strongly from low to high moments: For n = 1 the 
contributions from narrow resonances and continuum are 
roughly comparable, for n = 3 the continuum contribu- 
tion amounts to about 10%. Furthermore, the experi- 
mental contribution to the error decreases with increas- 
ing n from 9 MeV to 5 MeV, the /^-dependence, reflecting 
the theory uncertainly, increases from 2 MeV to 7 MeV. 
Despite the significant differences in the composition of 
the errors, the results are perfectly consistent. Since the 
result from n = 1 has the smallest dependence on the 
strong coupling and the smallest total error we take as 
our final value 



m c (3 GeV) = 986(13) MeV, 



(11) 



and consider its consistency with n — 2, 3 and 4 as ad- 
ditional confirmation. 
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n 


m c (3 GeV) 


exp 


OLs 


A* 


np 


total 


1 


986 


9 


9 


2 


1 


13 


2 


976 


6 


14 


5 





16 


3 


978 


5 


15 


7 


2 


17 


4 


1004 


3 


9 


31 


7 


33 



TABLE III: Results for m c (3 GeV) in MeV obtained from 
Eq. The errors are from experiment, a s , variation of \i 
and the gluon condensate. 



n 


m 6 (10 GeV) 


exp 


a s 


/' 


total 


m b (m b ) 


1 


3597 


14 


7 


2 


16 


4151 


2 


3610 


10 


12 


3 


16 


4163 


3 


3619 


8 


14 


6 


18 


4172 


4 


3631 


6 


15 


20 


26 


4183 



TABLE IV: Results for m 6 (10 GeV) and m b (m b ) in MeV 
obtained from Eq. Q. The errors are from experiment, a s 
and the variation of fi. 



Transforming this to the scale-invariant mass 
Tn c (m c ) [25|, including the four-loop coefficients of 
the renormalization group functions one finds [29I ] 
m c (m c ) = 1279(13) MeV. Let us recall at this point that 
a recent lattice determination, combining a lattice simu- 
lation for the data for the pseudoscalar correlator with 
the perturbative three- and four-loop result 11. 17. 2rj] 
has led to m c (3 GeV) = 986(10) MeV [13] in remarkable 
agreement with [l[ and the present analysis. 



The same approach is also applicable for the case of 
the bottom quark. Using the new moments with their 
significantly reduced experimental error (see Tab. HU, one 
obtains the results for the bottom quark mass at the scale 
\i = 10 GeV as listed in Tab. IIVI In comparison with the 
previous determination a minute upwards shift of 1 MeV 
(resulting from an upward shift of +3 MeV from the new 
data and a downward shift of —2 MeV from the new 
theory input) and a reduction of both experimental and 
theory error is observed. The three results based on n = 
1, 2 and 3 are of comparable precision. The relative size 
of the contribution from the continuum above 11.24 GeV 
which is modelled by perturbative QCD decreases for the 
higher moments n — 2 and 3. On the other hand the 
theory uncertainty, exemplified by the [i dependence is 
still acceptable. We therefore adopt the result from n = 2 
(which is roughly between the n = 1 and n = 3 values 
and which exhibits the smallest error) as our hnal result 



m fe (10 GeV) 
m b (m b ) 



: 3610(16) MeV, 
: 4163(16) MeV. 



(12) 



These values are well consistent with the previous deter- 
mination [H m 6 (10 GeV) = 3609(25) MeV and m h {m b ) = 
4164(25) MeV. 



It is straightforward to evolve the new value for m b to 
the normalization point at Mz and mt(mt) = 161.8 GeV 

m b {M z ) =2835 ± 13 ± 17 MeV , 
m 6 (161.8 GeV) =2703 ± 12 ± 19 MeV, (13) 

where a matching to the nj = 6 flavour theory has been 
performed in order to arrive at m(,(161.8 GeV). The first 
error originates from Eq. (fl~2|) the second from Sa s . 

For some of the applications it might be useful to ex- 
plicitely exhibit the a s dependence of our result, which 
is given by 



m c (3 GeV) 
mi, (10 GeV) 



m b (m b ) 



m b (M z ) 
TO fc (161.8 GeV) 



(986- 
(3610- 
(4163- 
( 2835- 
( 2703- 



0.1189 



0.002 
a, - 0.1189 



• 9 ± 10 MeV 





0.002 




- 0.1189 




0.002 


a s 


- 0.1189 




0.002 


u s 


- 0.1189 


0.002 



12 ± 11 MeV 



12 ± 11 MeV 



27 ± 8 MeV , 



2s : s ) MeV , 
(14) 



where a s — a s (Mz)- When considering the ratio of 
charm and bottom quark masses, part of the a s and of 
the [i dependence cancels 



m c (3 GeV) 
m b (10 GeV) 



0.1189 



0.002 



0.2732 — 0.0014 ±0.0028, 

(15) 

which might be a useful input in ongoing analyses of bot- 
tom decays. 



V. SUMMARY 

Based on new four-loop results for the higher deriva- 
tives of the vacuum polarization function and new 
BABAR data for bottom quark production in the thresh- 
old region, a reanalysis of the charm- and bottom-quark 
mass determination has been performed. The new data, 
a posteriori, give additional sup por t to the analysis of 
CLEO data presented in Ref. [2CJ and, furthermore, 
lead to a significant reduction of the experimental er- 
ror. The new theory results for the higher moments 
lead to a further reduction of the theory uncertainty and, 
equally important, demonstrate the consistency between 
the analysis based on different moments. The final re- 
sults, m c (3GeV) = 0.986(13) GeV and m 6 (10GeV) = 
3.610(16) GeV are consistent with the earlier determina- 
tion in Ref. [l| and, together with Ref. j^, constitute the 
most precise determination of charm- and bottom-quark 
masses to date. 
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